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~ ■ Abstract 

M-t' 
1) ■ 
*—, In this work we revisit Wald's cosmic no-hair theorem [1] in the context of accelerating 

Bianchi cosmologies for a generic cosmic fluid with non-vanishing anisotropic stress 

tensor and when the fluid energy momentum tensor is of the form of a cosmological 

constant term plus a piece which does not respect strong or dominant energy conditions. 

Such a fluid is the one appearing in inflationary models. We show that for such a system 

anisotropy may grow, in contrast to the cosmic no-hair conjecture. In particular, for 

a generic inflationary model we show that there is an upper bound on the growth of 

o 

the upper bound is found to be of the order of slow-roll parameters. We examine 
our general discussions and our extension of Wald's theorem for three classes of slow- 
roll inflationary models, generic multi-scalar field driven models, anisotropic models 
IT/i ■ involving U(l) gauge fields and the gauge-flation scenario. 
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1 Introduction 

The Universe at cosmological scales and as we see it today looks homogenous and isotropic. 
To understand this almost perfectly homogenous and isotropic Universe, given that the 
cosmic evolution may start from a generic initial condition over which we have no control, 
it is natural to seek a "dynamical" explanation: isotropic and homogenous Universe is an 
attractor of the cosmic evolution. According to the standard model of cosmology cosmic 
evolution is governed by the Einstein gravity coupled to the cosmic fluid, this latter idea if 
true, should be an outcome of Einstein equations. The first such attempt was made in [2j [3] 
arguing that the late-time behavior of any accelerating Universe is an isotropic Universe. 
This statement was dubbed as "cosmic no-hair conjecture" . The first attempt to prove this 
conjecture was presented in Wald's seminal paper pp. Wald's cosmic no-hair theorem states 
that Bianchi-type models (except Bianchi IX) with the total energy momentum tensor of 
the form T^ v = —Aog^ + T^, and with a constant positive A and a T^ v satisfying Strong 
and Dominant Energy Conditions, respectively SEC and DEC, will approach de Sitter space 
exponentially fast, within a few Hubble times if -1 = a/3/A . For a more thorough historical 
review and other related works on cosmic no-hair conjecture/theorem see [4] and references 
therein. 



Cosmological observations, specially the CMB data [5], indicate that the early Universe 
have in fact gone through such an accelerated expansion inflationary period. As a result, if 
cosmic no-hair conjecture/theorem holds, all traces of initial anisotropy and inhomogeneity 
should be washed away [6j. On the other hand, inflation has ended and we are not in a 
de Sitter Universe with Hubble parameter equal the one at the end of inflation, as Wald's 
theorem suggests. This is simply because inflationary setups are constructed to be so and 
do not strictly respect SEC or DEC assumptions of Wald's theorem. (For a more rigorous 
discussion of the latter fact see appendix B.) Cosmic no-hair theorem should hence be re- 
visited for inflationary models. Such studies has of course been carried out extensively, e.g. 
see [7] and references in jl] for an incomplete list. It may still seem plausible to expect that 
the result of Wald's theorem, namely an isotropic Universe, to hold for inflationary models. 
This is due to the fact that inflation is expected to have lasted around 60 e-folds, i.e. 60 
Hubble times, while the time scale for anisotropy damping, if Wald's theorem holds, is a few 
Hubble timeso 

The above expectation, and applicability of cosmic no-hair conjecture for inflationary 
models, may be challenged in two different ways: by changing the setup upon which Wald's 
theorem is based. That is, modifying gravity theory used for inflationary model building; 
or alternatively, to search for particular inflationary models which drastically violate SEC 
or DEC assumptions. Both of these ways have been studied, e.g. see [8] for the former and 
[9j [10] for the latter. All of our analysis in this work will be within Einstein gravity and 
hence applicable to the latter class of models. 

As pointed out and will be discussed briefly in section 3 and in more details in appendix 
B, during inflation assumptions of Wald's theorem are violated. However, this violation 
should be such that the model allows for a long enough inflationary period and has viable 
observational signatures. It has been shown that it is indeed possible to construct such stable 
anisotopic inflationary models, e.g. see [TUl HI [12] • There is hence a need to revisit the fate 
of Wald's theorem in more general inflationary setups and in particular slow-roll models. 
This is the question we tackle in this work. 

In this work, in section 2, we study evolution of Bianchi cosmological models within 
Einstein gravity sourced by cosmic fluid with anisotropic stress. As in Wald's work p] we 
decompose the energy momentum tensor as — A(t)^ All/ + 7^ v , with A > and T Ml/ satisfying 
SEC and WEC, but now allow for A to have time dependence. As proved in appendix 
A, generic inflationary models conform to this decomposition. In section 3, we push our 
analysis further by concentrating on inflationary frameworks, and in particular slow-roll, 
quasi-de Sitter expansion. This allows us to refine our results, obtain an upper bound and 
study the time evolution of the anisotropy. In section 4, we apply our general arguments of 
sections 2 and 3 to two classes of models: scalar-driven inflationary models and inflationary 
models with background gauge fields. In the last section we summarize our results into 
the "inflationary extended cosmic no-hair theorem" and briefly discuss the observational 
consequences of our results. 



1 As we will discuss in sections 3 and 4, this argument may be used for cases where the energy momentum 
of the system driving inflation does not contribute to anisotropic stress, like scalar-driven model. 



2 Extended cosmic no-hair, the general setup 

The Bianchi family are cosmological models with spatially homogeneous surfaces which are 
invariant under the action of a three dimensional Lie group, called the symmetry group. A 
Bianchi model can always be written as (e.g. see |13j ) 

ds 2 = -dt 2 + hije* <g> e j , (2.1) 

where i,j = 1,2,3, label the coordinates in homogeneous space-like hypersurfaces S t and 
hij = hij{t) is the spatial metric which can be decomposed as 

hij = e 2a e 2 ^. (2.2) 

Here e a is the isotropic scale factor and /3y is a traceless matrix which describes the anisotropics. 
Furthermore, e l are characterizing one-forms with the following property 

de* = -^cV^Ae fc , (2.3) 

where c % - k are the structure constants of the corresponding Bianchi type. 

Let n^ be the unit tangent vector field of the congruence of timelike geodesies orthogonal 
to the homogeneous space-like hypersurfaces S t . Then, we obtain the following covariant 
form for the spatial metric h^ u (12. ip 

hp, = g^u + n^n u , (2.4) 

where g^ is the metric of the space-time. The extrinsic curvature of S t is defined as 

•*»■ fJ,V = T\ n V- v = ~n IM 1 \ ) 

where the dot represents derivative with respect to the proper time t. One can decompose 
K^ u into trace, and trace- free parts 

K^ = -Khfu, + a^ , (2.6) 

where 

K{t) = K^hT = 3H(t) . (2.7) 

Here H(t) = a is the Hubble parameter corresponding to the homogenous scale factor e a . 
The shear of the time-like geodesic congruences a^ is related to the time derivative of /3y 
and is a symmetric, traceless and purely spatial tensor: 

h, u a^ = n^ = 0. (2.8) 

We now analyze Einstein's equation for Bianchi models 

G>„ = T^ = -A{t) 9llu + V, (2.9) 



where A(t) is a positive but time-dependent cosmo logical term and in our conventions we set 
8nG = 1. We assume that 7^ u satisfies the strong and weak energy conditions. It is shown 
in the appendix A that it is always possible to decompose energy momentum tensor of any 
inflationary system in this way. The strong energy condition (SEC) states 

(T^ - -7g^)tH u > 0, for all time-like f. (2.10) 

The dominant energy condition (DEC) stipulates that 

7 fill t fJ, t' v > for all future-directed causal vectors t^,t' u . (2.11) 

We note that the above for t ,fl = t M leads to weak energy condition (WEC) c J llu t^t u > 0. 
One can then decompose Einstein equations (12. 9p into four constraint equations 



V«"^ = ->R--a^ + -K 2 -A(t), (2.12) 

T^V = K\c\^K%c\ x , (2.13) 

and six dynamical equations 

{%x ~ l -h <TX 7)h%h x j = Z n (Jii + \{Kh tj + K 2 h i:j + Ka i:j ) - 2a iX a) - \{t)h tj + ( %- . (2.14) 

(3) 

Here Rij is the spatial Ricci tensor and can be written in terms of the structure-constant 



tensor c* fc (12. 3p as 

'X = 7 C ^ C / - c Vfe) - c ^ c l > ( 2 - 15 ) 



(3) D - 1 r, r, kl „k „ I „ Akl) 



where indices raised and lowered with metric hy, and R = R^h^ is the spatial curvature 
of E t . Note that all but one of Bianchi models have negative spatial curvature 

(3) i?<0; (2.16) 

Bianchi type IX has a positive curvature R > [TJ [13]. 

Combining trace of (12. 14j) with (I2.12p . we obtain Raychaudhuri equation 

(V - \g^)n»n v = -K + A(£) - ^K 2 - a^ v , (2.17) 

where K = £ n K. Contracting (I2.14p with W and removing the trace part, we have the 
following equation for the shear tensor 



a\ + Ka% + } S% = 7 kl h ki h\ - -7 k ih kl ti , . (2.18) 



where S % • is the anisotropic part of the spatial 3-curvature 



( V i = ( V i - ~ (3) i?/iV (2.19) 



Integrating equation (12.181) . we obtain the following integral equation for a 1 At) 

a i. {t) = e -Mt) f ^21 (T uh ki h i - hruhWj - %) + Cy^ , (2.20) 

with integration constants C J -. In order to determine a 1 At) we need more information about 
the energy momentum tensor 7^ u . 



3 Extended cosmic no-hair theorem for inflationary backgrounds 

Equations of previous section was written for quite general A(t) and 7^, and although we 
restricted our discussions to A(£) > and to 7^ v satisfying SEC and DEC, these conditions 
were not used. In this section, we focus on our main interest, inflationary systems and 
crucially use these conditions and investigate the dynamics of shear tensor and the evolution 
of anisotropics during inflationa 

Inflation is defined as an epoch in the history of universe in which the scale factor has 
an accelerating expansion 

(e a ) 

±-^- = H + H 2 = H 2 (l - e) > 0, (3.1) 

where e(t) = — -§% is a positive and growing quantity and inflation eventually ends when 
e > 1. Although condition ( 13. ip concerns the dynamics of isotropic part of the geometry, we 
argue here that inflation strongly restricts dynamics of anisotropic part of metric and puts an 
upper bound on the growth of anisotropies. Note that in spite of some similarities between 
our set up and pQ, the behavior of these two systems are totally different. The increasing e(t) 
which is a generic property of inflationary models, is in contrast with the dynamics described 
by the Wald's theorem [I], in which e(t) is a positive quantify that asymptotically approaches 
zero. More technical details about the differences between these two set ups might be found 
in Appendix B. 

Without loss of generality 7^ may be decomposed as 

V(t) = p(£K^ + P(t) V + *>(*), (3-2) 

where n^ is the anisotropic stress tensor which is symmetric, traceless, and purely spatial: 

n ia ,h ta, = & 7^71" = 0. 

We also define A = 3Hq, where subscript denotes the initial value (at the beginning of 
inflation) and 

A(t) = A + 5A(t) . (3.3) 



2 In the cosmology literature, it is more common to work with the Hubble parameter H(t) instead of 
K(t) (which is more common in GR literature). For the rest of this work we will use H(t) which from (|2.7[) 
is equal to \K. 



From the combination of (j2.12p and (12.171) . after inserting (13.21) and (13. 3p . we obtain 

m-P(t) +5m = ^ r + h-A + 3H \ (3.4) 

m±m = b-R-H- 1 -^. 0.5) 

2 6 2 M y J 

Since H is a decreasing quantity during inflation, from ( 13. 4 p we learn that, ^ +<5A(t) is a 
negative quantity in inflationary models of all Bianchi types, expect type IX. Nonetheless, the 
spatial curvature R has a time dependence proportional to e~ 2a and this term is damped 
quickly in inflationary systems. That is, R is damped after a few e-folds and becomes 
negligible. Therefore, for all Bianchi models, including Bianchi type IX, even if 2 + 
SA(t) is not negative to start with, it becomes negative very quickly. After a few e-folds R 
is negligible and we approximately have 

5A(t) p{t)-P(t) _ H(t) , A , _ e(t) | H 2 



3# 2 (£) 6# 2 (£) 3H 2 (t) 3H 2 (t) 3 # 2 (£) 

c(t). (3.7) 



p(t)+P(t) , cr^cr^ ff(*) 



2# 2 (t) 2# 2 (£) # 2 (£) 

On the other hand, T Ml/ satisfies SEC which choosing t^ = n^+s^, where s M is a normalized 
arbitrary space- like 4- vector normal to n M , implies 

p(t) + P(«)+7r^(t)sV>0, «,j = 1,2,3, (3.8) 

where s l = s^K 1 ^. Recalling the fact that -k^ is traceless and the above inequality should be 
hold for all s\ one obtainqj 

KWI<_4^ = 4 £W V,, (3.9) 



H 2 {t) ~ H 2 (t) 

From the combination of (13. 6p and (13. 7p and noting that (13.81) implies p(t) + P(t) > we 
obtain 

a^a^ < 2H 2 (t)e(t) , 
M W <-(p W + ^)<0. <3 - 10) 

The above already implies a weak upper bound on anisotropy 

[cr%-I <H(t)V2e. (3.11) 

As we will see in the next section assuming slow-roll leads to a stronger bound. 

During inflation (on the average) e is an increasing quantity and this result among other 
things shows the possibility of growth for the anisotropic part of stress tensor. This is in 



'Although we have not used here, we remind that WEC on T M „ imply p > 0. 



contrast with behavior of systems described by the Wald's theorem, in which all elements of 
T^, e.g. 7T%(t)'s, are damped exponentially with a time scale Hq 1 pp. 

Inserting (j3.2p into ( I2.20p . we have the following form for the Hubble-normalized shear 
tensor, ^ 

H(t) ~ 3H(t) J to de H(t>) + C > e ' (3 ' 12) 

where C 1 ^ is a constant at most of the order y/e ( 13. lip and S l j(t) oc 0(— H)e~ 2a ^>. There- 
fore, the second and the third terms in (13.121) have a damping behavior and soon become 
negligible. Then, after a few e-folds, we have the following relation for the Hubble-normalized 
elements of shear tensor 

**,•(*) ^ 1 [" d{H{t')e^')) vrV(t') 



#(i) 3H(t)e^J io (1 + J^) # 2 (i')' 

where t is a few e-folds after the beginning of inflation. As we see g | ffl acts as a source 
term and its five degrees of freedom may be determined for specific models of inflation. Some 
examples will be reviewed in section 4. 

For models with perfect fluid T^ U) in which 7T%(£) is identically zero, shear tensor a 1 At) 
is damped exponentially fast. Thus, in this kind of systems inflation washes away any initial 
anisotropies and isotropizes the system, in accord with the cosmic no-hair conjuncture. This 
is not necessarily the case in a general inflationary model with non-zero anisotropic stress 
tensor 7i l j(t). 

To summarize, we have so far studied the dynamics of anisotropies during inflation in a 
model independent way. We showed that in presence of anisotropic stress tensor n 1 ■ elements 
of shear tensor a 1 At) are allowed to grow by inflationary dynamics, which is against the 
cosmic no-hair conjuncture. Although a 1 At) can increase, as one expects from (13.91) . inflation 
enforces an upper bound value on their enlargements. 



3.1 quasi de Sitter expansion 

Thus far we considered a general model of inflation. In the following we focus on the slow-roll 
models. We investigate the dynamics of system assuming that n l j(t)/H 2 (t) always saturates 
its maximum value, and determine the possible upper bound value of anisotropies during 
slow-roll inflation. 



We start with defining the so called slow-roll parameters □ 



-£• " = -^' (314) 



which during the slow- roll are both very small and almost constant. 



Note that our definition of slow-roll parameters e and r\ and those usually used for the single scalar FLRW 

jip 2 - V(tp), e = i^y and 1)4, = ^, 



inflationary theory with £ = ^ip 2 — V(ip), e^ = |-w and ry^ = K^, are related as e^ = e and Tjj, = r) + e 



From (I3.5P and (13.91) . we find that 7^ u , a^a^ and R are at most of the order eH 2 . 
Hence recalling (1B.3J) and using f 1 2 . X 2 [) . we obtain 

H{t) ~ H (l-e H t), (3.15) 

where ~ means to the first order in slow-roll parameters, subscript denotes the initial value 
and H = a/A /3. Inserting the above relation in (13. 6p and (13. 7ft . we obtain u 

^ +£ l^r = - e ° (1+ww) ' (316) 

^? + i^y - £o ( 1 + 2(£o -, oW ). (3,7) 

As we see in (13.171) . o^ v o^ lv is not directly related to the slow-roll dynamics, but the combi- 
nation of 2H 2 <t) + 2H2(t) should be slow varying during inflation. As a result, a^ is not 
fully determined by the slow-roll dynamics, while (I3.16P implies that the isotropic part of 
the system is governed by the slow-roll. Thus, assuming slow-roll inflation, anisotropies can 
still evolve quickly in time. 

Applying slow-roll approximation in (I3.13p . we obtain the following form for the late time 
behavior of the Hubble-normalized shear tensor ,/,,, 



Hit) 



H (l + e H t)e-^ [ dt>e 3Hot '(l - 2e H tf)^J-. (3.18) 

J to H (t) 



In order to determine J, ■. , we need 7r J (t)'s five degrees of freedom which should be pro- 
vided by independent equations. However, regardless of the particular inflationary model 
which we consider, recalling (13. 9p . slow-roll enforces an upper bound on the enlargement of 
anisotropies. Rewriting (13. 9 j) in terms of slow-roll parameters, we obtain 



\^j(t)\ 
H\t) 



< Ae{t) ~ 4e (l + 2(e - 770)^0*)) • (3.19) 



As mentioned before, from the combination of (13.71) and (13.13!) . we realize that during the 
slow- roll inflation, ir % At) / H 2 it) can saturate its upper bound and grow in time. Then, 
assuming that ir l j(t)/H 2 (t) always saturates its upper bound during the slow- roll, we can 
find an upper bound value on the enlargement of a 1 jit) /Hit) enforced by the slow- roll 
dynamics. 

Inserting (13.191) into (13.181) . we obtain the upper bound value on the enlargement of 
anisotropies at the end of slow-roll inflation 

l(Tjl \t sl <l(eo-Vo)- (3-20) 



H ' sl ~ 3 



Note that we are considering simple slow-roll models where e(£) is an always increasing function during 
inflation, e > 0. From the definition of e (I3.14J) we learn that -4- = 2(e — 77). 



where t s i is the end of slow- roll inflation and approximately is Ho€ot s i — 1. Thus, a gen- 
eral model of slow-roll inflation by the end of slow-roll inflation, leaves the cosmos almost 
isotropic, with a shear of the order eo- Note that considering slow- roll has improved the 
bound 1ETTT1) to (1X201) . 



4 Anisotropy dynamics in two classes of inflationary models 

The discussions of previous sections was made for generic inflationary models. In this section, 
we consider some examples and discuss the possibility of saturation of the anisotropy upper 
bound found by specific classes of slow-roll models. In particular we study scalar-driven 
models and models involving vector gauge fields in the inflationary background. 

4.1 Scalar driven inflationary models 

In this part with in the context of Bianchi cosmology, we investigate the late time behavior of 
anisotropy in three scalar models of inflation: ordinary multi-scalar field model, K-inflation 
and DBI inflation. 

• Ordinary multi-scalar field models with the action [6] 

S = J rf 4 xv^(| - \d^jd^j - V{fa)), (4.1) 

where I runs from 1 to N and (pi's are N scalar fields, 

• K-inflation with the action [B] 

S = J d 4 x^=g(* + P(<p, VvO) , (4.2) 

where tp is a scalar field and X := |(V</?) 2 , 



DBI inflation with the action [14 



S = j d'x^-g [| - j (V© - 1) - Vtf] 



(4.3) 



where D := det(<5^ + fGijd IJ 'ip I d v ip J ), I runs from 1 to N, ipi's are scalar fields and 
Gu is an internal metric which determine the interaction between yj/s. 

The characteristic of all the above scalar driven models is that their energy momentum tensor 
T^y is of the form of a perfect fluid and hence the shear stress 7t(t) l j is identically zero in these 
systems. Equation ( 13. 12ft then implies that in all of the above multi-scalar driven inflationary 
models the anisotropy a^ u damps out exponentially fast in a couple of Hubble times, leaving 
an isotropic background. These models respect the inflationary no-hair conjuncture. 

9 



Note that due to the correspondence between f(R) gravity models and general relativity 
with a scalar field matter, we expect these systems to show similar behavior in anisotropy 
damping (regardless of the details of f(R)) and hence to respect the inflationary no-hair 
conjecture. 

4.2 Models of Inflation involving vector gauge fields 

In this part we consider two classes of models with a vector gauge field turned on in the 
background level. As far as anisotropies are concerned, these models are more interesting 
than scalar-driven cases. In both of the models, as we will see, the energy momentum tensor 
has an anisotropic stress n % • which can source anisotropies. Despite this fact, the dynamics 
of the two models is such that anisotropy does not grow in one class of models. 

To illustrate a detailed analysis for these models we restrict ourselves to Bianchi type 
I model, while the results seem to be generic to all Bianchi models. Bianchi type-I axially 
symmetric metric can be described by the line element 

ds * = _ dt 2 + ( t )2( e -M*) dx 2 + e Mt) {dy 2 + dz 2)) ; (44) 

where a(t) is the isotropic scale factor and e a ® represents the anisotropy. Due to the vectorial 
nature in these models, energy-momentum tensor has a non-zero anisotropic stress 7i(t) 

r; = diag{-p{t), P(t) - Tr(t), P(t) + ^vr(t), P(t) + ^(t)), (4.5) 

which its dynamics (whether it increases or decreases during inflation) determines the fate 
of anisotropies during inflation. Note that from ( 14.4)) . the shear tensor a^ u and a are related 
as 

a^ v = 6a 2 , (4.6) 

hence here instead of working with a^ v , we use a. 

4.2.1 Inflationary universe with anisotropic hair 

In [10j[T5], M. Watanabe, S. Kanno and J. Soda introduced and studied an inflationary model 
with anisotropic hair, as a counterexample to the cosmic no-hair conjecture. Motivated from 
supergravity, their model includes scaler field(s) as inflaton field(s) coupled to a massless 
U(i) gauge field. The action of the model is given as 



S — I d A X\f—g 



^-la^V-yfoO-I 



- -d^d'ip - V{y) - -J 2 {^)F, U F v 



(4.7) 



where /(</>) is the coupling function of the inflaton field to the vector one and the field 
strength of the vector field is given as F^ u = d^A u — d u A^. Below we review the analysis of 

urn Ha. 



10 



Choosing the temporal gauge, A = 0, homogeneous fields are taken to be of the form 
v4 M = (0, A x (t), 0, 0) and tp = tp(t). With this ansatz, one can solve equation of motion for 
the vector field as 

A x = r 2 Me- 4 -% (4.8) 

a(t) 



where Pa is an integration constant. Substituting the ansatz and C I4.8[) into the action, we 
obtain the energy density p(t), pressure density P(t) and anisotropic stress n(t) 

Note that for this specific model the anisotropic stress 7r(£) is equal to \pa, where pa is the 
contribution of the vector field in the energy density. 

Let us consider the chaotic inflation with the following potential 

V{<p) = ~mV, (4.9) 

and choose the coupling function to be f(ip) = e cip I 2 . After inserting V((p) and f((p) 
into (14. 9p . one can solve Einstein equations. Being interested in cases in which the energy 
density of the (gauge) vector field grows during inflation, c > 1, we realize that the Hubble- 
normalized anisotropic stress (J^- = — ) increases in time. As a result, (13. 12[) enforces the 
Hubble-normalized shear j| to grow during inflation. 

In [10], the evolution of the Hubble-normalized shear a/H for cases with c > 1 has 
been calculated numerically. They found two slow-roll phases and as expected, the Hubble- 
normalized shear grows during inflation. In FigJTJleft panel, we present their results for a/H 
which is plotted for various values of parameter c under the initial condition \fcipi = 17. As 
wee see, in the first slow-roll phase, the solutions show a rapid growth. During the second 
phase it still grows but with a slower growth rate. In the second slow-roll phase, in which 
a ^. Ha, we approximately have 

m^im ( 4io) 

H ~ 3 H* ' { ' 

i.e. the Hubble-normalized shear grows as a result of the growth of Hubble-normalized 
anisotropic stress during the inflation. As has been discussed in [TTJ , it is possible to choose 
the initial conditions and parameters such that the model has only one phase during about 
60 e-folds of inflation. This, as shown in the right panel of FigfTl happens for fi = UMpi 
and c = 2. As we see, the Hubble normalized anisotropy grows during inflation and gradually 
saturates its upper bound value and becomes of the order e. 

These models exhibit an anisotropic attractor with non-zero, nevertheless small a/H. As 
depicted in Fig.l, at the end of inflation a/H ~ e, saturating our bound ( 13.20)) . We comment 

11 
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Figure 1: In the left panel we have the evolutions of Hubble- normalized shear 
various c values when \f~cxpi = 17. As we see, growing during inflation and becoming of 
the order e, j| saturates the upper-bound value. This figure is taken from [10]. Then, in 
the right panel and for a different initial values and parameters, we have the evolution of 
Hubble-normalized shear jj '•= jj during inflation when c = 2 and = ipt — UMpi. Again, 
anisotropies grow during inflation and eventually saturate their upper bound value. This 
figure is taken from [TT] . 



that although the two-phase model of [10] does not strictly satisfy our slow-roll conditions 
for the two phases, since the second phase lasts long enough (c/. Fig. 1), this model obeys 
our general arguments and the upper bound. 



The above analysis has been extended for more general Bianchi models in [12J and for 
cases with more scalar fields and non-Abelian gauge fields [TSJ US]. In all these cases the 
anisotropy dynamics is in accord with our general analysis of section 3, and can saturate the 
anisotropy bound (13.201) toward the end of inflation. 



4.2.2 Non-Abelian gauge field inflation, gauge-flation 

In [IT], we introduced a novel inflationary scenario, non-Abelian gauge field inflation or 
Gauge-flation for short. In this model inflation is driven by non-Abelian gauge field min- 
imally coupled to Einstein gravity. It was shown that non-Abelian gauge field theory can 
provide the setting for constructing an isotropic and homogeneous inflationary background. 
This was achieved noting that gauge fields are defined up to gauge transformations and 
that any non-Abelian gauge group has an SU(2) subgroup. The global part of this SU(2) 
subgroup can be consistently identified with the rotation group. In particular, this means 
that among the gauge field components A ^ (// being the space-time and a the gauge index) 
we have turned on a specific combination which behaves as an scalar under rotation. This 
scalar combination is coupled to all the other components of the gauge field through gauge 
interactions and hence one may excite these other components either classically or quan- 



12 



turn mechanically, violating the isotropy and rotational symmetry of the background FLRW 
trajectory. 

As discussed in details in [17], one particularly convenient choice of the gauge- fiat ion 

action is 

~R 1 „„ „ ..„ K 



S = Id Xy/—g 



fT a p ft" _i ( fivXarpa rpa \2 

a [iv a ' oo A \ \M> Act/ 



(4-11) 



where ^ vXu is the totally antisymmetric tensor and the gauge field strength tensor F a is 
given by 

F%, = d,A\ - d v A\ - ge\ c A\A c v) (4.12) 

here fi,u,... labeled the space-time indices and run from to 3, while a,b,... labeled the 
indices of the SU(2) gauge symmetry algebra and run from 1 to 3. Our action consists of 
a Yang-Mills term (which due to scaling invariance can not lead to an inflating system) as 
well as a specific F 4 term, Tr(F A F) 2 , which its contribution to the energy- momentum 
tensor has the equation of state p 4 = —P 4 , perfectly suited for driving an almost de Sitter 
expansionjj 

Stability of the FLRW inflationary trajectory in gauge-flation against small fluctuations 
around the isotropic FLRW background was already studied and established in [T7J. Then, 
in [T5] which will be briefly reviewed here, starting from Bianchi type-I cosmology and 
through analytical and numerical studies, it was shown that the isotropic FLRW inflation is 
an attractor of the dynamics and the anisotropics are damped within a few e-folds, in accord 
with the cosmic no-hair conjecture. To do so, let us start with fixing the temporal gauge 
A a = 0, with the ansatz 

A = A a i T a dx i = e -^) ^w) T i & + e a{t) X{t)a{t)^{t){T 2 dy + T 3 dz) , (4.13) 

A (t) 

where T a are the SU{2) gauge group generators, [T a ,7V] = ie c ab T c . Note that A 2 (t) = 
1, a(t) — corresponds to the isotropic background of [T7] . 

It turns out that the equations take a simpler form once written in terms of 

<P(t) = a(t)i>(t) . (4.14) 

Substituting the ansatz and the axi-symmetric Bianchi metric into T^ v for the gauge fields, 
we obtain p(t), P(t) and ir(t) as 



P = Pk + P YM , P = ~ Pk + -Pr MJ 



1 

(4.15) 



6 As discussed in [17] and in more details in |18j . from particle physics model building viewpoint, a 
Tr(F A F) 2 type term can be argued for, considering axions in a non-Abelian gauge theory and recalling the 
axion-gauge field interaction term Laxion ~ ^Tr(F A F) 2 . Then, integrating out the massive axion field ip 
leads to an action of the form we have considered. 
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where p K and p YM are respectively contributions of F A and Yang-Mills terms 
3 Kg 2 <f) 4 2 



Pym 



2 
3 
2 



a 4 a 

1 
3A 4 



' (^ Of , A N<^ 2 , 2 ^ 2 /0 . /. * 0x2 . 

( 2 (or + -)-) + _(- + (<7-~)~) + 



c/ 



A a' 



A a' 



(2 + A 6 ) g 
3A 2 



2^4 



(4.16) 

(4.17) 



and 



a 



(A" 4 (A e 



2 ) 



2a 2 (3 + A- 4 (2 + A 6 )g) 

Note that p K is only a function of <fi(t) and not A(t). As we see, in the isotropic case A 2 = 1 
(a = and A = 0), ir(t) vanishes. 

Assuming that the system undergoes a quasi-de Sitter slow-roll inflation, the dynamics 
of A and a was examined and shown, both analytically and numerically, that this system has 

two attractor solutions A = {^f) 3 = ±1 which regardless of the initial values of A, all the 
solutions converge to them within a few e- folds |19j . These two attractor branches, which 
are physically identical due to parity and charge conjugation invariance of our gauge-nation 
action (14. lip , correspond to the isotropic quasi-de Sitter solutions. Thus, gauge-nation is 
globally stable with respect to the initial anisotropies and respects cosmic no-hair conjecture. 



H 1 
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Figure 2: Hubble- normalized anisotropic stress tensor, -^J-, and Hubble- normalized shear 
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M£Z := £W f or a sys tem with k = 3.77 x 10 15 , g = 10" 1 , ip = 0.6 x lO" 3 , ip = 10 
A = 10, A = —3.6. As we see ^f increases (exponentially) for a very short time in the 
early stage of inflation (in small H t), saturating our upper bound (I3.20p . Then, it is damped 
exponentially fast to its isotropic fixed point. 



The slow-roll isotropic inflation is the attractor solution in our system, nonetheless, due 
to the vectorial nature in gauge-flation, for very large and small initial values of A there is 
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the possibility for ■£ to saturate its upper bound value (13.201) for a very short lapse of time in 
early stages of inflation. Both our numerical and analytical calculations reveals that in the 
extreme limits of A 6 3> 1 and A 6 <^C 1, there is a region where anisotropy grows exponentially 
for a very short period, before getting exponentially damped to its isotropic fixed point. 
Although in these cases gauge-nation does not follow the strict dynamics indicated by the 
cosmic no-hair theorem jl] for the very short period of time, the anisotropies are indeed 
damped within few Hubble times, in accord with the cosmic no-hair conjecture. 

Although the analysis of [19] was carried out for Bianchi type-I model, a similar behavior 
for all Bianchi models is expected, once we assume having a quasi de Sitter slow-roll expansion 
phase. This is due to the fact that this inflationary phase is driven by p K which dominates 
the p YM and that the anisotropic stress ir(t) only receives contributions from the Yang-Mills 
term (TCT5]) . fl416|) and fHTTj) . 

5 Concluding remarks 

In this work we extended Wald's cosmic no-hair theorem for general inflation setups, and in 
particular for slow-roll models. We proved an "inflationary extended cosmic no-hair theorem" 
which states that despite the fact that anisotropies can grow in time during inflation, their 
amplitude remains small and there is an upper bound on their amplitude. 

For Bianchi cosmological models and in the context of general relativity, we investigated 
the dynamics of anisotropies during inflation. Considering a general imperfect fluid form 
for the energy-momentum tensor T^ u we showed that the late time (after a few couple of 
e-folds) behavior of anisotropic tensor a^ v is governed by the anisotropic stress tensor 7^; if 
Tti{t) is vanishing or damped in time, as ( I3.12p shows, the shear follows the same behavior 
as 7r* •. Choosing the orthonormal frame in which T^ v is diagonal, a^ u may have off-diagonal 
elements in general. However, since these elements have no source term they are quickly 
damped, with time scale Hq. 

Here we proved that inflation enforces an upper bound on Hubble-normalized anisotropic 
stress ir l j(t)/H 2 (t), proportional to the slow-roll parameter e(t), which can in principle grow 
in time during (slow- roll) inflation. As a result, inflationary dynamics allows the dimension- 
less anisotropies <j\(t)/ H(t) to grow during the inflation, in contrast to the cosmic no-hair 
conjuncture. Furthermore, assuming slow- roll inflation, the dynamics of anisotropy does not 
necessarily follow the slow-roll dynamics and they can generally evolve quickly. This effect 
was first noted in [10] in the context of a special inflationary model with Bianchi type I met- 
ric in general relativity. Although, inflation allows for the growth of the diagonal elements 
of dimensionless anisotropy, as we showed here it puts an upper bound \a l j\/H < |(e — r] ) 
( 13.201) at the end of slow-roll inflation. 

Our work was motivated in part by possible observational prospects and traces that a 
small but non-zero primordial anisotropy may have on the CMB, which may show up as 
statistical anisotropy in the CMB power spectrum. Currently, there are only bounds on the 
statistical anisotropy, for some studies in this direction see e.g. [15l HEl 120] . These bounds 
will hopefully be improved by the upcoming cosmological observations and may be used as 
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a tool to pin down on the model of inflation. 

Acknowledgement 

We would like to thank Jiro Soda for his collaboration in the early stages of this work 
and for comments on the draft. The work of AM is supported in part by the fund from 
Boniad-e Mell Nokhbegan of Iran. 

A Useful decomposition for energy momentum tensor of inflation- 
ary systems 

Here, we prove that it is always possible to describe the energy momentum tensor T^ v of any 
inflationary system (which satisfies DEC, but violates SEC) as 

T^ = -A(t)^ + V, (A- 1 ) 

where A(t) > 0, and 7^ u which satisfies both SEC and WEC. 
Assuming inflationary dynamics in a system, we have 

H + H 2 > , 

which using ( 12. 17p gives the following inequality for the total energy-momentum tensor T^ v 

(T M „ - \g,jr)n»n v + o^ = -3(H + H 2 ) < 0. (A.2) 

As a necessary condition for the above inequality to be satisfied, T^ u should violate SEClH 
However, decomposing T M „ as (1A.1I) . one can choose A(t) in such a way that 

A(t)>3(H + H 2 ) + a^, (A.3) 

and obtain 

(V - l -g^)n»n» > 0. (A.4) 



That is, choosing A(t) as (1A.3J) one can obtain a 7^ v which respects SEC. 



On the other hand, for 7^ v to satisfy WEC, we need 7 IIU n fl n' / > 0, which from (j2.12p 
requires 

A(t)<3H 2 + ^R-±a^. (A.5) 

Then, demanding 7^ v to satisfy both SEC and WEC, ( 1A.3I) and (1A.5I) imply that A(t) should 
be in the following region 

3(H + H 2 ) + cv<r^ < A(t) < 3H 2 + ^R - ^a^. (A.6) 



We comment that having inflation in generic anisotropic models demands a condition stronger than 
violation of SEC, it requires (T M „ — ^g IMV T)n ,J "n 1 ' < —cr^cr^, as indicated in (|A.2[) . 
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In order to prove the validity of such a choice, one should show that the upper bound value 
for A(t) is less than its lower bound value. In the following by showing that the difference 
of the upper and lower bound values is a positive quantity we prove the validity of ( 1A.6I) . 



The total energy-momentum tensor T^ v satisfies DEC which implies that T fiU n fl n u + 
\T^W U > 0. We hence have 

which recalling ( 12. 12ft and (" 12. 1TI) yields 

-3H-~a^ v + ^R>0. (A.7) 

The above is exactly the difference of the upper and lower values in (1A.6J) . We have thus 
proved the validity of the statement made in the beginning of this appendix. 

We note that as is implicit in the above, SEC and WEC conditions does not uniquely 
specify A(t). One may hence assume extra conditions on 7^ v . One such choice is to take 
T M „ to be stiff matter, that is a matter with p = P (with possibly non-zero anisotropic stress 
71^). This freedom and various choices of A(t) does not lead to a stronger bound on the 
anisotropy. 



B Wald's theorem assumptions and inflationary models 

Consider inflationary models in which we have exponential expansion in a finite period of time 
and inflation eventually ends. These systems do not respect conditions of Wald's theorem 
[Tj. Systems which are described by Wald's theorem evolve toward the de Sitter solution 
and inflation will never stops in them. Here we review in more technical details exactly how 
models of inflation fail fulfilling assumptions of Wald's cosmic no-hair theorem. 

To this end, we employ Wald's notations p] where total energy momentum T^ v (12 .9p was 
decomposed in terms of a positive cosmological constant Ao, and an extra term T^ v 

T^ u = -Aog^ + T M „ , (B.l) 

with T satisfying DEC and SEC. This latter is not respected in inflationary models. In terms 
of our notations, 

f, w = -(A(t) - A )<7 M „ + V • (B.2) 

As we find from the above equation for any given T^, there is an ambiguity in defining T^ U) 
unless the value of A is also specified. In order to resolve this ambiguity, we choose A to 
be 

Ao = 3# 2 L , (B.3) 



which using (12.1 2D . determines initial value of n^n^Tf 



I 1 " 



n"n% v \^ o = i ( X, - IcvoH^, (B-4) 
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where So represents spatial hypersurface of initial time at the beginning of inflation. Recall- 
ing that for all Bianchi models expect type IX, R < 0, equation (1B.4J) implies that T f 



M^lEn 



does not respect the dominant energy condition in these models. 

Noting that H 2 is a decreasing quantity in inflationary systems, and using (12.X2H . ( IB. 21 ) 
and (IB. 41) . one finds out that 

• for all Bianchi cosmological models expect type IX, we have n^n^T^ < during the 
inflation, which implies that T^ v violates DEC. 

(3) 

• In addition, in case of having inflation in Bianchi type IX in which R > 0, the spatial 
curvature is rapidly damped and after a few e- folds, we have n^n u T IMU < 0, which again 
implies T^ u violates DEC. 

• Since during inflation H 2 is a decreasing quantity, one can show that it is not possible 
to find a A > in such a way that the resulting T^ u satisfies SEC and DEC for the 
whole period of inflation. 

Note that although T^ eventually violates DEC in inflationary systems, the total energy- 
momentum tensor T^, always satisfies DEC and the energy density of the system, T^n^n", 
is a positive definite quantity. To summarize, we showed that during inflation, if T^ u satisfies 
SEC it will necessarily violate DEC, which is against the requirements of Wald's cosmic 
no-hair theorem. 
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